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A general formalism to solve nonlinear differential equations is given. Solutions 
are found and reduced to those of second order nonlinear differential equations in 
one variable. The approach is uniformized in the geometry and solves generic nonlin- 
ear systems. Further properties characterized by the topology and geometry of the 
associated manifolds may define global properties of the solutions. 



1 Introduction 



The systematic solution to non-linear partial differential equations has prohibited 
many advances in mathematics and physics. These equations, contrary to standard 
theory and linear equations, appear disparate and unsolvable in the general case. In 
the past the solutions to these systems of equations have involved many equation 
dependent techniques that are different in various regimes. A systematic formalism 
to the solutions of these problems is required. 

In this presentation a foundation is given based on geometric methods allowing 
for a uniform treatment of the solutions to these sets of non-linear partial differential 
systems. The treatment here is to give a solution set to these equations and to provide 
an optimum calculus for future improvements. 

The general non-linear system is solved for by finding the geodesies on an asso- 
ciated manifold. Even in the generic multi-dimensional system, the solution space is 
found by solving only second order non-linear equations in one proper time variable 
(e.g. the geodesies). The reduction in dimensionality is useful for computations and 
for manifesting global properties of the solutions, such as existence and behavior in 
various regimes. 

2 Setup 

The work begins with the algebraic solution to systems of polynomial equations as 
generated in [T], involving toric and geodesic properties as found for example in [21, [Sj, 
[3]. The sets of equations are those of the form, 



with c labeling the equation, and <J c (j, I) labeling the Ith term of the cth equation. 
The permutations of the terms labeled by a. An example set is, 



rn 





z n x + 2z™zl + 3z£ = 



(2.2) 



3z 2 m + 2z%z™ + z? = . 



(2.3) 
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The equations in (|2.1|) are general and contain, for example, the well known case 
pertaining to Fermat's last theorem, 

P(zi) = z? + z$-z5 (2.4) 

Equations with m = oo are also contained in the set. 

The individual terms are accorded to the function(s) solved for together with 
the derivatives. For example, in the Fermat's set of equations, one would have Z\ = 
u(w,t), z 2 = u x (w,t), and z 3 = u t (w,t); this pertains to a quantum mechanics 
problem in the two-dimensional space of (w,t). 

The interest here, however, is in modeling general sets of (non-linear) differential 
equations in a uniform manner. In the interpretation of (|2.1|) . the dynamics of the 
system is governed by a set of equations with initial conditions spanning a (potentially 
singular) manifold. To set a basic groundwork, the equations in (j2.1)l parameterize a 
holomorphically complex manifold as in pQ. 

The spaces considered are (toric) Calabi-Yau, although more general manifolds 
are required. These manifolds have varieties defined by the equation set in (j2.1|) . and 
may potentially be singular. The beginning setup is then, to every non-linear differ- 
ential equation, there is an equation describing a fibration of a Calabi-Yau manifold 
Mz over the possibly curved spacetime Mx- The fibration is general and an example 
is, 

P(z; x)=z™ + R{x)z™ + 4 = 0, (2.5) 

with R(x) a generic function describing the x dependence of the non-linear differential 
equation. The other coordinates are z\ = u(x,t), z 2 = u x (x,t), and z 3 = u t (x,t). 
Hence, the fibration is obvious in this example, as is the differential equation. Most 
general forms are described by, 

p c {z i ;x j ) = Y,a*(j,) n Qm( x ) z *m • ( 2 - 6 ) 

j=l 1=1 

These equations are modeled by the flow in a 'product' of Z eg) X, and are analyzed 
in this work. 

One example in particular is the Navier-Stokes equations, 
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which has a manifold interpretation in terms of Vi = Z(i,o), v (i,j) = If 1 with j = 
0, 1, . . ., and p(x) and fi general functions in the space X. The remaining equations 
is contained in the divergence condition, 



£fr =a (2 - 8) 

These equations are solvable in the approach here. 

The fibration Z <g> X is important to consider first. First, the initial conditions 
are defined by a point in this space. 

The solutions are found via flows in the total space X®Z . There are integrability 
conditions in the direction of the flow given the derivatives, i.e. u(x) = Z and 
u x (x,t) = d x Z. As an example, consider a manifold Mz spanned by the variables, 



which are functions of the base space variables x and t. Then the dynamics is de- 
scribed by the differential equation P(u,u x ,u xx ,Uxxt), with u(x,t) a function of the 
variables on the base R(x,t) (e.g. flat space R 2 ). Label the functions as, 



du dzi dz 2 dz 3 

Z!=U, z 2 = — = z 3 = — — , z A = — . (2.10) 

X ox ox ot 

The initial conditions u°, u°, u° xx1 and u° xxt at (x°, t°) generate a point on the fibration 

Mz <8> X . (These initial values must satisfy the differential equations.) 

In order to find the values of u (and u x , u xx , u xxt ) at a different point (x,t), the 
former variables follow a flow from the point (x°, t°) to (x, t) while preserving the 
conditions in (|2.10j) . The path in Mz generated through the equations in (|2.10|) are 
integrability conditions on the coordinates (x°,t°) to (x,t). 



3 The case of P(z t ) = 0. 

There are two cases of interest, when P(zi]Xj) = P{zi) and when there is explicit base 
dependence. The former case is analyzed first. The Zi coordinates are the functions u, 



u XJ and etc. The values from the initial point z° to the final point Zj are constrained 
on the manifold P(zi) and are generated via a geodesic flow. Then the coordinates Xj 
are determined via the integrability conditions in (j2.10j) ; the inclusion of an auxiliary 
complex space V{zi) is generally required. 

The first step, as in pQ, is to model a space via the polynomial equation P(zi) = 
(with possibly the non-holomorphic modification of V(z). The geodesies on these 
space will require solutions. 

Formally the spaces used in the work are toric Calabi-Yau, which may have singu- 
larities; the general manifolds are holomorphically toric, that is, have two polynomials 
specifying the both the holomorphic and nonholomorphic sides. A modified M = 2 
D-term specification can be used to find these metrics. 

These polynomials give information about the existence of solutions, further from 
the locus of points in (|2.1jl . In the Kahler examples, label the space pertinent to 
(J2.1J) as Mp( Zi ) and its Riemannian metric as g^. Its Kahler so that both g^ u = 
d^d u \Q.(f){zi,Zi) (in terms of z and z, g = g^) and T p ^ u = Xjldpd^dy In <p(zi, Z{) hold. 
The polynomial sets of finite degree are modeled by a finite dimensional Calabi-Yau; 
the equation sets of infinite degree (i.e. transcendental) are described by an infinite 
dimensional manifold. 

The geodesic equation is, with the coordinates x = (zi,Zi), 



d 2xP _|_ pp.^i/ dXfj, dx v _ q (3 1' 

dr 2 dr dr 

The Kahler example admits a complex form with the connection components, 



l/2^^1n0(x M ,^) . (3.2) 



The coordinates x contain both the holomorphic and anti-holomorphic pieces describ- 
ing the geometry. Its complex Kahler form is 



didz^ + dp_g] d i H(l)) \ _ ^Igjgi^ = o 
dr dr dr dr 2 



(3.3) 



or 




Figure 1: The solution via flow in the space Mz <8> Mx- 

These equations are second order in derivatives, but non-linear because of the Kahler 
potential. Solving these equations gives the z\ coordinates as a function of the 'proper 
time' r; the 2n initial and final coordinates z° and z{ determine the unknowns. 

Once the path from z° to Zi is found, the coordinates in the base are determined. 
As the Zi coordinates are explicit functions of the initial conditions and r, the in- 
tegrability conditions in (j2.10j) are first order in derivatives and the coordinates Xj 
requires solutions in r. Take the example of the equation, 

u 5 + u\ x + u x = . (3.5) 

The Zi are Z\ = u(x), Z2 = u x (x) and z% = u xx (x). The integrability conditions in 
(j2.10j) requires, 
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The function x is found via, 



dx = 1 dzy , s 

dr z 2 {r) dr ' 1 ' 



The coordinate in ()3.7|) can be solved by an integral 



dx = -^dr (3.8) 

and gives its r dependence. Further equations need to be solved, and these put 
conditions on the solutions to Zj. For example, the equation z 3 = dz 2 /dx is consistent 
only if, 



Zl ' T - Z2 - T _^ 1 d _ /o 
zi z 3 2 or or 

The equation in (J3.9j) is satisfied only if the final coordinates z(, and z( (and the 
initial z°, z%, and z 3 at the point x) are chosen well. These further requirements are 
algebraic in nature, and are hand in hand with the existence of solutions. 



4 Integrability conditions 

In determining the coordinates Xj from the conditions in (j2.10j) there are complica- 
tions due to the number of equations and the number of free variables X{. These 
complications are removed by enlarging the toric space to include a non-holomorphic 
side P(zi) = 0, which may be of a large degree. The manifold in principle is no 
longer Kahler, and holomorphically separated. The solution to the non-linear pde's 
is solved by the geodesic flow represented in the holomorphic coordinates Zj, while 
the constraints are satisfied via the influence on Z{ via the flow in the non-holomorpic 
side Zi. This is analyzed in this section. 

Consider a differential equation containing, 

Z\ — U Z2 = U X1 Z 3 = U X2 Z4 = U XlXl Z 5 = U X2X2 Zq = u XlX2 . (4.1) 
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There are fewer coordinates Xi than integrability conditions in (|4.1|) . After solving 
for the coordinates x\(t) and X2{t) with, 

dxi " dx2 

there are four integrability conditions, 

dz2 dz2 dz3 dz% , . 

3 — — ^4 -i— = z s -r- = z 6 — = z 5 , (4.3) 
ctxi dx2 ax i dx2 

and two relations between them, 

dz2 dz 3 dz 2 dzj, 

Za Zq Zq Z$ 

The integrability conditions, e.g., 



(4.4) 



d/Z2 y u3?i 

dr dr 



z 4 (4.5) 



determine four constants, i.e. four of the z{ in solving the geodesic motion. However, 
if the solutions are not possible, which is very probable, then the final points are 
taken to be arbitrary and a different route (or geodesic) is required; this different 
route changes the functional form of the solution Z{ so as to make these integrabibility 
conditions satisfied. The change in the route is made possible by altering the manifold 
on the non-holomorphic side P(z), while preserving the holomorphic polynomial P{zi) 
which models the differential equation. 

In a separate issue, perturbations may be added to make the integrability condi- 
tions linear in derivatives; this makes the analysis easier in solving these equations. If 
the differential equation does not contain the first derivative in the coordinate, then 
either: (1) the first non-vanishing term (e.g. u xx or higher order) must be used to 
define the spatial coordinate Xj, or (2) a perturbation of the form u x has to be added 
to the differential to define Xj as in the previous. The first option requires a second 
order differential equation to solved, as opposed to a first order one, z 2 = d 2 zi/dx 2 , 



, . o d cc . dec , o / _ \ 

* = ^ + ( a7» (4 ' 6) 

Equations of these types (to be solved for x) are more complicated, and first order 
linearity is preferred. The perturbations otiU Xi are simpler to add, and post the 
analysis the coefficients aij are taken to zero. Of course, if the solutions are singular 
then the previous method should be used. 
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The example given in ()3.5|) describes how the general non-linear differential equa- 
tion is solved for Zi(r;z°,Xj) and Xj(r; z°,x°). The general solution requires two 
non-trivial steps: (1) the solution of a non-linear second-order differential equation 
describing the geodesic condition, and (2) the solution to a set of algebraic conditions 
describing the integrability relations. These steps are uniform to all of the non-linear 
partial differential equations that satisfy P(zi) = 0, i.e. without any Xi dependence. 

The integrability conditions, if not satisfied via the solution to the direct equations 
of motion via the toric variety with complex conjugates P{zj) = 0, and P{zj) = 0, 
then a procedure is developed based on a different manifold is used. Consider a 
general space defined by the constraint, 

p c (^)=E^«n^ (i) > ( 4 - 7 ) 

<x(i) i 

and 

£a*W II . ( 4 -8) 

a(i) i 

with possible logarithmic modifications ln(z) and ln(z). In general, the Laurent ex- 
pansion with positive and negative powers of the coordinates span all of the non- 
analytic terms. There are in principle an infinite number of constants in the general 
non-holormpohic side represented in (|4.8|) . These constants are used to satisfy the 
integrability conditions. 

The point of introducing the coefficients a a and a a is to allow functional depen- 
dence in the components of the Christoffel connection, so as to allow the integrability 
conditions (e.g. (|2.10jl ) to be satisfied. In general the metric is a function of these 
coefficients a and a. The holomorphic piece P(zi) = is used to describe the an- 
alyzed non-linear differential equation. The non-holomorphic constraint P(zi), with 
the coefficients a have enough terms (an infinite defining a non-analytic function) so 
as to allow the constraints to be satisfied. 

As the holomorphic side of the manifold is parameterized by P(zj), the final 
solution for = Z^irf) will be obtained by a geodesic condition in the total manifold 
spanned by P(z{) and P(z~i, in such a manner that the consistency conditions are 
satisfied. The latter are maintained by choosing the P(zi) and the initial and final 
non-holomorphic coordinates, z^ and Zf, appropriately to repeat their solution via the 
geodesic motion in the total space. The solution to the consistency conditions could be 
complicated due to the possibly complicated solutions to the coordinates Zi, although 
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in only one variable. The geodesic equations are non-holomorphic in the sense that 
both Zi and zf are required (and the mixed indexed T) in the individual flows; 
however, the solution to the non-linear partial differential equations are obtained 
from only the holomorphic piece. 

In the procedure described, the general metric on a non-complex (holomorphically 
toric) manifold is required together with the solution to the geodesies on the space. 
However, the holomorphic side is described by the original function P{zj) describing 
the non-linear partial differential equation. The reproduction of the integrability 
conditions, given for the one example in ([2. 10)) . requires algebraic solutions to the 
parameters (the solutions may require transcendental algebra). 

The primary difficulty in the solution is solving for the geodesies in the Calabi- 
Yau manifold and also solving the algebraic equations. However, the treatment of all 
the non-linearities is reduced to only second-order partial differential equations in one 
variable (an analysis of algebraic equations is presented in [T], requiring also second 
order partial differential equation solutions). 

The existence of solutions and the singularities are described by the steps in 
the previous example, and are generally given by the behavior of the metric on the 
associated manifold. It is relevant that the non-linearities and the types of solutions 
may be described by the geodesies and properties of the underlying geometry. The 
consistency conditions describe a submanifold connected to the initial conditions z° 
and x°. 

The solutions via this method generate the r dependence in the coordinates Xj 
and the functions z i: i.e. Xj(r; x°, zfy and Zi(r; x°, z%). The r dependence may be 
solved for as a function of Xj and substituted in the solution to find u (u x , . . . ), or 
the Zi, in order to manifest the solution's Xj dependence. 

The interpretation of the auxiliary non-holomorphic components Zj may be put in 
the form of integrability constants, perhaps with the the non-holomorphic polynomial 
P(z) having a group theory interpration. 

Spatial dependence of the initial conditions should be commented on. The initial 
conditions of u(x, t), u x (x, £),..., correspond to a point in the manifold parameterized 
by P c (zi). A spatially dependent 'wave-packet' at time t = as an initial condition 
may be constructed by choosing an appropriate point Q in the manifold and finding 
the values in a localized region of space x G Mi at time t = in which u(x, 0) is given, 
via the transports from the point Q to the 'final' points in this region. The solutions 
at these 'final' points are evolved into the required regions t ^ and the x. 
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5 The case of P[z{\ xj) = 0. 



The primary difference in this case is the inclusion of explicit coordinate dependence 
in the set of coupled non-linear partial differential equations. For example, there 
might be a general term p(x) as in the Navier-Stokes equations, 

|-+x;^-A".-^+/.(m). (5.D 

More general equations could have the x coordinates spanning a curved spacetime, 
with the (coupled set of) non-linear equations taking values on it and there being 
explicit metric dependence. 

The general x dependence in the equations described by P a {zi] Xj) is incorporated 
by enlarging the flow in the Zi coordinates to include Xj. In this manner, the con- 
straints P a (zi]Xj) are always preserved. The non-holomorphic polynomial(s) P a (z t ) 
describe the remaining side of the manifold and are used in the geodesic flows as in 
the previous case of no explicit x-dependence. 

6 Concluding remarks 

The general set of coupled non-linear differential equations is described and reduced 
in terms of flows in manifolds, and the solution of algebraic equations. The solution 
to the differential equations is described by geodesic flow, that is they reduce to 
second order non-linear differential equations in one variable only. The multi- variable 
aspect is further transferred to the required solving of a set of algebraic equations (of 
potentially infinite degree) in an auxiliary set of variables a. 

The second order different iality in one variable bears a resemblance to a La- 
grangian particle description. A Hamiltonian first order description would be useful. 

The nonlinear partial differential equations describes a manifold; due to the type 
of manifold, the solutions inherit modular properties and potentially further ones 
akin to a mirror symmetry on toric varieties. Furthermore, the manifolds' topological 
properties potentially may be used to classify solutions. 

The geodesic flows, and the sub-manifolds they parameterize, describe chaotic 
versus periodic solutions and characterize the solutions according to initial condi- 
tions. For example, a closed geodesic obviously generates non-chaotic behavior; these 
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topological properties characterize the integrability, and ergodicity, of the solutions. 
An ideal situation entails cohomological calculations to show integrability properties. 

The obvious use of the proper-time solutions to the non-linear partial differential 
systems, and their form, merits further investigation. Also, the connection to the 
geometry of the manifolds perhaps requires input from its topological properties for 
characterizations of the solutions. The explicit form of the geodesies, and the metrics 
required in order to find them, are not presented in this work. 
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